Conductance distributions in disordered quantum spin-Hall systems 
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We study numerically the charge conductance distributions of disordered quantum spin-Hall 
(QSH) systems using a quantum network model. We have found that the conductance distribution 
at the metal-QSH insulator transition is clearly different from that at the metal-ordinary insulator 
transition. Thus the critical conductance distribution is sensitive not only to the boundary condition 
but also to the presence of edge states in the adjacent insulating phase. We have also calculated 
the point-contact conductance. Even when the two-terminal conductance is approximately quan- 
tized, we find large fluctuations in the point-contact conductance. Furthermore, we have found a 
semi-circular relation between the average of the point-contact conductance and its fluctuation. 

PACS numbers: 73.20.Fz, 73.20.Jc, 73.23.Ra, 73.43.Nq 



I. INTRODUCTION 

The integer quantum Hall (IQH) effect is observed in 
two-dimensional (2D) electron systems in a strong per- 
pendicular magnetic field. 1-4 Both the two-terminal con- 
ductance and Hall conductances in the IQH insulating 
phase are exactly quantized in integer multiples of e 2 /h. 
Transitions between quantized values occur as the Fermi 
energy passes through the center of a Landau band. 
According to the theory of Anderson localization, all 
electrons in 2D systems with broken time-reversal sym- 
metry (unitary universality class) should be localized. 5 
Later, the quantization was understood to be related to 
the peculiar nature of the conducting states (chiral edge 
states), which are characterized by a topologically invari- 
ant Chern number. 6 The quantum Hall insulator is thus 
called a topological insulator to be distinguished from an 
ordinary insulator. 

Recently, the counterpart of the IQH effect in the sym- 
plectic symmetry class, the so-called quantum spin-Hall 
(QSH) effect, has been observed in 2D systems with 
strong spin-orbit interactions in zero magnetic field. 7-11 
The QSH insulator exhibits both quantized charge and 
spin-Hall conductances and is a topological insulator that 
is characterized by the Z2 topological number. 8 The pur- 
pose of this paper is to investigate the universal proper- 
ties (the properties that are independent of microscopic 
details) of systems which exhibit the QSH effect. To this 
end, we simulate numerically a suitable quantum network 
model. 12-14 We focus on the charge conductance, which 
is easily accessible in experiments. 11 

QSH systems have two kinds of insulating phases: a 
QSH insulating phase with edge states, and an ordinary 
insulating phase without edge states. In addition, the 
QSH system exhibits a metallic phase sandwiched be- 
tween the insulating phases. By tuning the sample setup, 



two Anderson metal- insulator transitions are observed. 11 
Numerical study of the critical exponent that character- 
izes divergence of the localization length 12,15 and of bulk 
multifractality of wave functions 13 indicate that the crit- 
ical properties of QSH systems are the same as those 
of conventional 2D symplectic systems, which is a uni- 
versality class with time-reversal symmetry and broken 
spin-rotation symmetry. 5 

In this work, we investigate the probability distribution 
functions of the two-terminal charge conductance for dis- 
ordered QSH systems. At Anderson transitions, the dis- 
tribution functions are sensitive to the universality class 
and the geometry including boundary conditions 16,17 of 
the system, but are independent of the model and sys- 
tem size. We show that the conductance distribution 
in the metallic and ordinary insulating phases, and the 
critical conductance distributions at the metal-ordinary 
insulator (M-OI) transition of QSH systems behave in 
the same manner as those of the 2D symplectic sys- 
tems. On the other hand, in the QSH insulating phase, 
the conductance clearly reflects the property of the edge 
states. Remarkably, the critical conductance distribution 
at the metal-QSH insulator (M-QSHI) transition is com- 
pletely different from that at the M-OI transition. It is 
to be noted that the critical conductance distributions 
depend not only on the symmetry and geometry but also 
on the presence of edge states in the adjacent insulating 
phase. In addition, we investigate the point-contact con- 
ductance, which illustrates the local transport properties 
of disordered QSH systems. Even when the two-terminal 
conductance is well quantized, the point-contact conduc- 
tance shows the fluctuating nature of the edge states of 
the QSH insulator. This indicates that the conducting 
state, which is almost free from back scattering, is not 
always confined to the edge. We also demonstrate a pe- 
culiar relation between the average and fluctuation of the 



point-contact conductance. 

This paper is organized as follows. We first describe 
the quantum network model for QSH systems in Sec. II. 
We next study the two-terminal charge conductance in 
disordered QSH systems in Sec. III. We then study the 
point-contact conductance in Sec. IV. The final section 
is devoted to the discussions and concluding remarks. 
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II. MODEL 

To describe the QSH system, we use a Z2 quantum 
network model, 12 > 13 > 18 which is based on the Chalker- 
Coddington model 4,19 that describes IQH systems. The 
Chalker-Coddington model is a regular network of di- 
rected current paths. It consists of links and nodes, which 
correspond to equipotential lines and saddle points of 
a random potential in the 2D plane, respectively. The 
Z2 network model can be constructed by introducing 
the spin degree of freedom into the Chalker-Coddington 
model so that each link has a Kramers pair, i.e., a pair 
of current paths for up and down spins propagating in 
opposite directions. The scattering of the Kramers pairs 
at an s-type node at a position i illustrated in Fig. 1(a) 
is described by the 4x4 scattering matrix Sj, 
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where Cj a is current amplitude at the link j with the spin 
cr(=t,|). 9^ m is defined as 
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where 8\ is a random phase associated with the link j 
for the node i, which is distributed independently and 
uniformly between [0,27r). The matrix s in Eq. (1) is 
defined through 2x2 matrices Q and I2 as 
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where the parameters p and q (s [0, 1]) correspond to the 
chemical potential and the strength of the spin-orbit in- 
teractions that mix up and down spins, respectively. The 
property of s and the choice of the phases in Eqs. (1) - (4) 
guarantee the time-reversal symmetry. The nodes next 
to the s-type nodes are the s'-type [Fig. 1(b)]. The scat- 
tering matrix s' t is defined by 90° rotation of the matrix 
Si to make the model isotropic. 

Let us explain boundary conditions of the Z 2 network 
model. In this work, we study conductance for the QSH 
system on which periodic boundary conditions (PBC) or 



FIG. 1. (Color online) Schematics of the (a) s-type and (b) 
s'-type nodes described by the scattering matrices Si and s£, 
respectively. The arrow on upper (under) side of a link indi- 
cates the direction of currents for up (down) spin. 
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FIG. 2. (Color online) Schematics of the Z2 networks in a 
square geometry with RBC The left and right figures cor- 
respond to the network model at p — > and 1, respectively, 
while the middle figure represents the metallic phase at inter- 
mediate p. When p is close to or 1, all the bulk states are 
localized in plaquettes. In the QSH insulating phase (p — > 1), 
however, two Kramers pairs of conducting channels run along 
the edges, which correspond to the helical edge state. 



reflecting boundary conditions (RBC) are imposed in the 
transverse direction. (Hereafter, PBC and RBC mean 
the boundary conditions in the transverse direction to the 
current if the direction is not explicitly written.) When 
RBC are imposed, due to the time-reversal invariance, a 
Kramers pair is totally reflected without spin-mixing at 
a node i located on the boundary, 
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We assume the arrangement of the nodes as shown in 
Fig. 2, where the QSH insulating phase is expected to 
appear for p — > 1, while the ordinary insulating phase 
appears for p — > 0. 

The system size of the network is measured in terms 
of the number of links; hence, a square system with size 
L contains L 2 links. In two-terminal structures for con- 
ductance measurements, L is equivalent to the channel 
number of Kramers doublets, which is assumed to be 
even. 

The detailed phase diagram of the Z2 network model 
has been examined by the finite-size scaling analysis in 
a quasi-one dimensional geometry 12 and is obtained as 
shown in Fig. 3. Because of its symplectic symmetry, 
the system shows a metallic phase sandwiched between 
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FIG. 3. (Color online) Phase diagram (Ref. 12) of Z2 net- 
work model with RBC. Note that the phase boundaries be- 
tween metallic and insulating phases are the same for net- 
works with PBC. For PBC, the ordinary insulating phase 
replaces the QSH insulating phase. Circled letters indicate 
the points considered in the corresponding subsections in 
Sec. III. The actual values of parameters (p, q) are as fol- 
lows: A: (0.365,0.500), B: (0.039,0.146), B': (0.783,0.146), 
C: (0.188,0.309), D: (0.561,0.309), and D': (0.563,0.146). 



two insulating phases. When PBC are imposed on the 
Z2 network model, both insulating phases correspond to 
the ordinary insulating phase. On the other hand, when 
RBC are imposed, the insulating phase located for larger 
p exhibits edge states and becomes the QSH insulating 
phase as expected. Thus, the Z2 network model shows 
three types of transitions: the M-OI transition for PBC, 
the M-OI transition for RBC, and the M-QSHI transition 
for RBC. 



III. TWO-TERMINAL CONDUCTANCE 

In this section, we consider two-terminal charge con- 
ductance for the various phases in the QSH system with a 
square geometry (L links in both the x- and y-directions). 
When the ideal leads are attached to the left and right 
terminals, the charge transport is described by a 2L x 2L 
scattering matrix S, 
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where V^rT denotes the incoming (outgoing) current 
amplitude on the left (right) lead, and t and t' (r and 
r') are L x L transmission (reflection) matrices. By em- 
ploying the Landauer formula, 4 conductance G in units 
of e 2 /h is given by 

G = Tr(t f t). (7) 

We have obtained the t matrix by using the transfer 
matrix method, 4 in which the transfer matrices can be 
constructed from the scattering matrices Sj and s[ in 
Eqs. (1) - (3) along with boundary conditions. 




FIG. 4. (Color online) Conductance distributions P(G) in the 
metallic phase A with RBC, for L = 64, 128, and 256 (from 
left to right). Distributions keep Gaussian shapes (dotted 
line) for increasing system size. 



A. Metallic phase 

First, we focus on the metallic phase. It is known that 
in the metallic phase of systems with symplectic symme- 
try, conductance distributions P(G) follow a Gaussian 
function with constant variance. 20 In addition, the aver- 
aged conductance (G) increases with logarithm of system 
size with a universal coefficient as a consequence of the 
anti- localization, 5 



(G) =G + 7r- 1 lnL/Z, 



(8) 



where Go (S> 1) is the Boltzmann conductance and I 
denotes the mean free path. 

We have calculated the conductance in the metallic 
phase of the Z2 network model at the fixed parameters p 
and q corresponding to A in the phase diagram (Fig. 3). 
We have confirmed that the distributions are well fitted 
by Gaussian functions. For example, the distributions 
for systems with RBC are shown in Fig. 4. We have also 
confirmed that the averaged conductance obeys Eq. (8) 
irrespective of boundary conditions (Fig. 5). Therefore, 
the metallic phase in the QSH systems is qualitatively the 
same as in the ordinary symplectic systems. It is noted 
that, recently, the anti-localization behavior in the metal- 
lic phase of the QSH system is analytically confirmed 21 
by using a tight-binding model proposed in Ref. 7. 



B. Insulating phases 

Next, we focus on insulating phases (B and B' in 
Fig. 3). It is known, for the ordinary insulating phase, 
that the averaged conductance exponentially decays as 
the system size increases and becomes zero in the ther- 
modynamic limit. In addition, the distribution function 
of conductance becomes a log- normal function, i.e., the 
distribution function of the logarithm of conductance, 
P(lnG), follows a Gaussian function. 

We have calculated the conductance in the ordinary 
(at B in Fig. 3) and QSH (at B') insulating phases. Fig- 
ure 6 shows the averaged conductances for each insulating 
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FIG. 5. (Color online) Averaged conductances as a function of 
system size in the metallic phase A with PBC (o) and RBC 
(□). The solid lines represent the anti-localization formula 
Eq. (8). 
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FIG. 6. (Color online) Averaged conductance as a function of 
system size, for ordinary insulator B with PBC (o) and RBC 
(□), and QSH insulator B' (x). 



phase rapidly decrease and converge to different limiting 
values with increasing system size L. Indeed, these limit- 
ing values of conductance indicate that the QSH insulat- 
ing phase has two edge states. The distribution func- 
tions of the logarithm of the conductance P(ln G) for the 
ordinary insulator are well fitted by Gaussian functions 
as shown in Fig. 7. On the other hand, for the QSH insu- 
lating phase, P(G) is approximated by the delta function 
S(G—2), which is an evidence of the edge states. The de- 
tailed properties of the edge states in the QSH insulating 
phase are discussed in Sec. IV. 



C. Metal-ordinary insulator transition 

Figure 8 shows the conductance distribution P(G) at 
the M-OI transition point (C in Fig. 3) in a square ge- 
ometry with PBC and with RBC. For comparison, we 
also calculated conductance distributions at M-OI tran- 
sition of the symplectic class by using a tight-binding 
model [the SU(2) model 22 ] with PBC and fixed bound- 
ary conditions (FBC) and a non-chiral network model 
(the S2NC model 23 ) with PBC and RBC, which do not 
show the topologically-nontrivial phase. It is found that 
the distribution functions at the M-OI transition calcu- 
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FIG. 7. (Color online) Conductance distributions P(lnG) in 

ordinary insulating phase B with RBC for L — 256 ( ) 

and with PBC for L = 128 ( — ). For ordinary insulating 
phase, the distribution functions of InG are well fitted by 
Gaussian functions (dotted line) for any system size. For QSH 
insulating phase, distribution functions P(G) asymptotically 
approach to the delta function S(G — 2). 



lated in the Z2 network model coincide reasonably well 
with those in the other models belonging to the sym- 
plectic class. We have confirmed that the system size 
dependence of the distribution functions is negligible for 
sufficiently large L and these results are expected to be 
reproduced everywhere on the line of the same type of 
transition. In fact, the same distribution has also been 
obtained at the transition point on the other side of the 
phase boundary (D) for PBC. This detail independence 
of the critical conductance distributions implies that the 
distributions are universal and the Z2 network model, 
i.e. the QSH system, belongs to the same Wigner-Dyson 
symplectic class as the SU(2) model and the S2NC model. 
This is consistent with the calculation of the critical ex- 
ponent of the localization length and multifractality 12 ' 13 



D. Metal-quantum spin Hall insulator transition 

Finally, we focus on the conductance distributions 
P(G) at the M-QSHI transition points (D and D' in 
Fig. 3). As shown in Fig. 9, the conductance distribu- 
tions at these transitions are completely different from 
that of the M-OI transition with RBC in Fig. 8. We em- 
phasize that these distributions are insensitive to details 
such as system size L (» 1) or the parameter (f>, q) (along 
the same phase boundary), and hence are considered to 
be universal. The distribution at the M-QSHI transition 
has peak structure around G = 2, while they are rather 
broad at the M-OI transitions (see Fig. 8). The positions 
of peaks seem to reflect the helical edge states in insulat- 
ing phases. These results are quantitatively clarified in 
Table I, which shows the average (G) and fluctuation, 



sdv(G) = ^a7(G) = V(G 2 > - (G) 2 (9) 

of the conductance for each transition. 



Q. 



C5 

Q. 




0.6 



RBC (FBC) 




q j rf WBWS M tfMfrwWMffWW* 




FIG. 8. (Color online) Conductance distributions at the M- 
OI transition point C calculated on the Z2 network model 

( ), S2NC model (o), and SU(2) model (x) for PBC and 

RBC [FBC for the SU(2) model]. More than 10 6 independent 
random configurations have been realized to draw each curve. 



TABLE I. Averages and fluctuations of the critical conduc- 
tance in units of e 2 /h. 

Transition (G) sdv(G) 

M-OI (PBC) 1.42 0.60 

M-OI (RBC) 1.27 0.65 

M-QSHI 2.37 0.35 




FIG. 10. (Color online) Distribution functions of the largest 

P(2ti) ( ) and second-largest P(2ts) ( — ) transmission 

eigenvalues for M-OI transitions C with PBC and RBC. Dot- 
ted lines are the distribution functions of the conductance 
P(G) in Fig. 8 for G < 2 and for G > 2, both of which are 
normalized to be 1. The latter is shifted by —2 along the 
horizontal axis to be compared with P(2ts). 



Here we focus on the largest and the second largest eigen- 
values T\ and T3. We also note that the transmission 
eigenvalues and the conductance are related as 




g = E- 



(12) 



FIG. 9. (Color online) Conductance distributions at the M- 
QSHI transition with L = 256. The distributions at the dif- 
ferent points D ( ) and D' ( — ) almost coincide with each 

other. 10 samples are calculated for each critical point. 



The transmission eigenvalue statistics gives us a bet- 
ter understanding of the results mentioned above. The 
transmission eigenvalue is defined as the set of eigenval- 
ues of tH in Eq. (7), {t 1 ,t 2 ,--- ,t l }, where we order 
them according to 



T\ > T 2 > •• • > T L . 



(10) 



We note that all of the eigenvalues are doubly-degenerate. 
Thus, 



T\ = Tl > T3 = T4 • ■ • > TL-1 = TL- 



(11) 



Figures 10 and 11 show the probability distribution func- 
tions of the largest and the second largest transmission 
eigenvalues P(2n) and P(2t3), for each type of tran- 
sitions. Like the IQH transition, the largest eigenval- 
ues determine the bare shape of the conductance distri- 
butions P(G) for the M-OI transitions with PBC and 
RBC. 24 The second largest eigenvalues contribute to the 
tails of P(G) for G > 2. For the M-QSHI transition, 
the doubly-degenerate largest eigenvalues are narrowly 
distributed near n = 1, which shift P(G) by about two 
along the G-axis, and the second largest eigenvalues de- 
termine the shape of P(G). This implies that well con- 
ducting channels exist even at the M-QSHI transition, 
although there is no edge state. These channels, which 
dramatically change the shape of P(G) at the M-QSHI 
transition, might be related to the QSH edge states in 
the adjacent QSH insulating phase. 



IV. POINT-CONTACT CONDUCTANCE 
A. Point-contact conductance 

To investigate the nature of QSH edge states, we have 
calculated the point-contact conductance. This is the 
conductance between two small probes (such as scanning 
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FIG. 11. (Color online) Distribution lunctions ol the largest 

P(2ti) ( ) and second-largest P(2ts) ( — ) transmission 

eigenvalues for M-QSHI transition D. Dotted lines are the 
distribution functions of the conductance P(G) in Fig. 9 for 
G < 2 and for G > 2, both of which are normalized to be 
1. The latter is shifted by —2 along the horizontal axis to be 
compared with P(2t3). 
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FIG. 12. (Color online) Schematics of the Z2 networks with 
point contacts (a) on the edge, y p — 1, and (b) in the bulk 
region, y p — L/2. Each point contact is connected with a 
link. PBC are imposed on the longitudinal direction. 



tunneling microscope tips). For the network model, the 
point-contact conductance is calculated as the conduc- 
tance between two links. 25 ~ 27 We consider a cylindrical 
geometry (2L links in the x-direction and L links in the 
y-direction) and two point contacts separated by distance 
L as shown in Fig. 12. The point-contact conductance 
depends on the position of the contacts in the y-direction. 
We assume that both contacts are located at the same 
distance from one of the edges. That is, the contacts are 
attached at (x,y) = (0, y p ) and (L,y p ). 



B. Method 

To introduce point contacts into the network, we cut 
link k and link I where the contacts are attached. We 



then define incoming current amplitudes (c 
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and outgoing current amplitudes (cjS*, c^YSc^*,^*) on 
the corresponding links. The current amplitudes satisfy 
the equation 
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where the 4L 2 x 4L 2 scattering matrix S for all links 
of a network consists of the 4x4 scattering matri- 
ces Si and s'j in Eqs. (1) - (4) for a node. For 
given (cS., cjS , cX, cX), the remaining current amplitudes 
(ci t ,---,clf,c%f,---,c^ t ,c^\---,c 2 L^) are uniquely 
determined by the following set of 4L 2 simultaneous lin- 
ear equation with 4L 2 unknowns 
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As a consequence of the structure of these equations, 
there is a linear relationship between the incoming and 
outgoing current amplitudes 
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The point-contact conductance G pc is given by 

^pc = I 1 ' l^pc^pcji 

in units of e fh. 



(16) 




FIG. 13. (Color online) Current amplitude (|c;-|-| 2 + |c^| 2 ) of 
a Z 2 network with (p,q) = (0.750,0.146), for L = 80. The 
contacts are attached at the upper edge. 




FIG. 14. (Color online) Distribution functions of the point- 
contact conductance at the upper edge, y p = 1 ( ), near the 

edge, y p — 3 ( — ), and at a distance on the order of the local- 
ization length, j/ p = 19 (... ), with (p,q) = (0.750, 0.146), for 
L — 80. 10 samples are realized. Although the two-terminal 
conductance is quantized for corresponding parameters, large 
fluctuations appear for the point-contact conductance. 




FIG. 15. (Color online) (a) Averages and (b) fluctuations of 
edge (D:L = 80, +:L=160) and bulk (o:L = 80, X:£ = 160) 
point-contact conductances with q — 0.146 as functions of p. 
The lines are guide to the eyes. The plateau of fluctuation of 
the bulk conductance corresponds to the metallic phase. 



D. Relation between average and fluctuation 



C. Conductance near edge states 

First we focus on the local property of the edge state 
in the QSH insulating phase. Figure 13 shows the spa- 
tial distribution of current for a Z2 network in the QSH 
insulating phase. The current spreads on the order of 40 
links from the upper edge. This is consistent with the es- 
timate of the quasi-lD localization length (= 19.1 ± 0.2) 
with width 160, PBC, and (p,q) = (0.750,0.146). In 
this case, the point-contact conductance is not quan- 
tized and shows large fluctuations. Nevertheless, the two- 
terminal conductance remains well quantized as shown in 
Sec. IIIB. This implies a significant suppression of back 
scattering, once the system width exceeds the localiza- 
tion length. 

Figure 14 shows probability distribution functions of 
the point-contact conductance P(G pc ) at y p — 1,3,19 
in the QSH insulating phase. The distribution P(G pc ) 
for 2/ p = 1 has a peak at G pc = 2 with a rather broad 
distribution. In addition, the distributions for y p = 3—19 
have broad tails toward G pc = 2. This indicates that the 
width of the conducting path associated with a QSH edge 
state is on the order of the localization length. 



We then see the point-contact conductance as a func- 
tion of p. Figure 15 shows the average and fluctuation 
of the edge (y p = 1) and bulk (y p = L/2) point-contact 
conductance for q = 0.146. When the bulk conductance 
becomes small, two-terminal conductance is well quan- 
tized, since the small bulk conductance means the small 
inter edge state coupling, which causes the back scatter- 
ing. 

Unlike the two-terminal conductance, the point- 
contact conductance shows the system size dependence 
even at criticality and, then, it is difficult to find any uni- 
versal property in the point-contact conductance. How- 
ever, a universal rule appears when we plot the fluctu- 
ations as a function of the average of the point-contact 
conductance for various p (Fig. 16); the data points fall 
on a single curve. The curve is approximated by a semi- 
circle, although there are small intrinsic deviations. The 
relation is independent of the details of the model and 
the geometry (p, q, L, or y p ). Indeed, we have also con- 
firmed this semi-circler relation for the systems with RBC 
in both the x- and y-directions. 
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FIG. 16. (Color online) Fluctuations as a lunction of the 
averaged point-contact conductance for edge (□ : L = 80, + : 
£ = 160), 2/ p = l, and bulk (o:L = 80, x :L = 160), y p = L/2, 
with q — 0.146. This shows a specific relation similar to a 

semi-circle ( ) between the average and fluctuations We 

have confirmed that the same relation is satisfied for different 
geometry including boundary conditions, other contact point 
(j/ p = 3), or different spin-flip rate (q = 0.309). Statistical 
errors are smaller than symbols. 



V. DISCUSSION AND CONCLUDING 
REMARKS 

In this paper, we have studied the transport proper- 
ties of QSH systems by numerical calculations on the 
Z2 network model. We have shown that the conductance 
distributions at transition points are sensitive to the type 
of transition, i.e., the presence or absence of edge states 
in the adjacent insulating phase, and found the universal 
conductance distribution for the M-QSHI transition. We 
have also shown that the universal conductance distribu- 
tions for conventional symplectic systems with PBC and 
RBC are reproduced in the M-OI transition of the Z2 net- 
work model. This is consistent with Ref. 13, where the 
boundary multifractality has been shown to be sensitive 
to the existence of the edge states. 

In the QSH insulating phase, the point-contact conduc- 
tance fluctuates remarkably near the edge, while the two- 
terminal conductance is well quantized. We have found 
a universal relation between the average and fluctuation 
of the point-contact conductance. The similar relation 
was suggested for the two-terminal conductance in IQH 
systems. 28 Motivated by this, we have confirmed the re- 
lation for the point-contact conductance in IQH systems 
(Fig. 17) and found it is closer to a semi-circle. We have 
also calculated the conductance for diffusive transport in 



a single-channel quantum wire based on the Dorokhov- 
Mello-Pereyra-Kumar (DMPK) equation. 29 ~ 32 The curve 
(Fig. 17) agrees with our data in the small conductance 
region, but deviates from the semi-circle for large con- 
ductance, which can be achieved only by the edge state. 
The semi-circular relation, therefore, may be a charac- 
teristic of the topological insulators. Recently, much ex- 
perimental progress have been made in developing local 
probes 33 ' 34 for IQH systems. We expect that this relation 
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FIG. 17. (Color online) Fluctuations as a function of the av- 
eraged point-contact conductance in the Chalker-Coddington 
model (IQH system). Plots for both edge (□ : L = 80, 
+ : L = 160), y p = 1, and bulk (o : L = 80, x : L = 160), 

j/p = L/2, agree with the semi-circular relation ( ). The 

result for the single channel DMPK equation (Refs. 29 and 
30) ( — ) agrees with the semi-circle only for small averaged 
conductance. 



is experimentally accessible since this relation is indepen- 
dent of the parameters and the boundary conditions. 

The above conclusions are mostly based on the net- 
work model, the details of which are different from re- 
alistic samples. However, since the properties calculated 
here are expected to be universal (see for example, Fig. 
8, where excellent agreement of the tight-binding model 
with the network model is demonstrated), we expect our 
results can be verified by experiments. 
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